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We find an exact formula for the thermally averaged cross section times the relative velocity 
(cru r ei) with relativistic Maxwell-Boltzmann statistics. The formula is valid in the effective field 
theory approach when the masses of the annihilation products can be neglected compared with the 
dark matter mass and cut-off scale. The expansion at x = m/T 1 directly gives the nonrelativistic 
limit of (nUrei) which is usually used to compute the relic abundance for heavy particles that decouple 
when they are nonrelativistic. We compare this expansion with the one obtained by expanding the 
total cross section a(s) in powers of the nonrelativistic relative velocity v r . We show the correct 
invariant procedure that gives the nonrelativistic average (a nr v r ) n r coinciding with the large x 
expansion of {<7U r ei} in the comoving frame. We explicitly formulate flux, cross section, thermal 
average, collision integral of the Boltzmann equation in an invariant way using the true relativistic 
relative u re i, showing the uselessness of the Mpller velocity and further elucidating the conceptual 
and numerical inconsistencies related with its use. 


I. INTRODUCTION 

While there are compelling evidences in astrophysics 
and cosmology that most of the mass of the Universe is 
composed by a new form of non baryonic dark matter 
(DM), there is a lack of evidence for the existence of new 
physics at LHC and other particle physics experiments. 
On the theory side, many specific models with new par¬ 
ticles and interactions beyond the standard model have 
been proposed to account for DM. 

Under these circumstances where no clear indications 
in favour of a particular model are at our disposal, the 
phenomenology of DM as been studied in a model inde¬ 
pendent way using an effective field theory approach, see 
for example [1-23]. 

Measurements of the parameters of standard model of 
cosmology [24, 25] furnish the present day mass density of 
DM, the relic abundance, flh 2 ~ 0.11 with an uncertainty 
at the level of 1%. Any model that pretends to account 
for DM must reproduce this number, which, on the other 
hand, sets strong constraints on the free parameters of 
the model. 

When the DM particles are weakly interacting massive 
particles that decouple from the primordial plasma at a 
temperature when they are nonrelativistic, the relativis¬ 
tic averaged annihilation rate (crv le i) can be well approx¬ 
imated by taking the nonrelativistic average of the first 
two terms of the expansion of a in powers of the nonrel¬ 
ativistic relative velocity. With u re i we indicate the rel¬ 
ativistic relative velocity and with v r the nonrelativistic 
relative velocity , as defined in B. To describe collisions in 
a gas, and in particular in the primordial plasma, the ref¬ 
erence frame that matters is the comoving frame (COF) 
where the observer sees the gas at rest as a whole and 
the colliding particles have general velocities V\^ without 
any further specification of the kinematics. 

It is thus desirable to formulate cross sections and rates 
in a relativistic invariant way, such that all the formulas 
and nonrelativistic expansions are valid automatically in 


the COF. Obviously, invariant formulas give the same re¬ 
sults in the lab frame (LF), the frame where one massive 
particle is at rest, and in the center of mass frame (CMF) 
where the total momentum is zero. We will see that the 
key for the invariant formulation is v Te \. 

On the contrary, in DM literature [26] instead of u re i 
it is used the so-called Mpller velocity v, see B. That 
this is incorrect was already discussed in Ref. [27] but 
papers using v continue to appear. The problem with v, 
which is not the relative velocity, is its non invariant and 
nonphysical nature, for it can take values larger than c. 

In this paper we first find an exact formula for (crt! re i) 
as a function of x = m/T calculated with the relativistic 
Maxwell-Boltzmann statistics. The formula is valid in 
the effective field theory framework such that the masses 
of the annihilation products can be neglected compared 
with the DM and the cut-off scale. For concreteness we 
work with fermion DM. We find the thermal functions 
corresponding to various interactions and in particular 
those corresponding to s and p wave scattering in the 
nonrelativistic limit which is given by the expansion at 
x 1. This is done in Section II, and A contains some 
mathematical results needed for the derivation of the ex¬ 
act formula and its asymptotic expansions. 

Then, in Section III, we present the correct invariant 
method for obtaining the same expansion by expanding 
the total annihilation cross section a (s) in powers of v r . 

We then discuss in Section IV the problems with the 
use v, while the numerical impact on the relic abundance 
of some incorrect methods employed in literature is eval¬ 
uated in Section V. 

B is preparatory for the whole paper: we remind how 
relativistic flux, cross section, rate, collision term of the 
Boltzmann equation and thermal averaged rate can be 
defined in the invariant way in terms of i; re i showing the 
uselessness of the Mpller velocity. 
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Figure 1. s and t channel annihilation diagrams reducing to 
the effective vertex corresponding to the lagrangian Eq. (1). 


II. EXACT FORMULA FOR THE THERMAL 
AVERAGE IN THE EFFECTIVE APPROACH 

We consider a DM fermion field \ that couples to other 
fermion fields ip through an effective dimension -6 opera¬ 
tor of the type 

£a = -^(xr a x)(V?r 6 ^). (1) 

The DM particles can be of Dirac or Majorana nature and 
have mass m, while ip are the standard model fermions 
or new ones. Here A a ,b are dimensionless coupling asso¬ 
ciated with the interactions described by combination of 
Dirac matrices T 0j b. A is the energy scale below which 
the effective field theory is valid. In the exact theory A 
corresponds to the mass of a heavy scalar or vector boson 
mediator that appears in the propagators. The ip masses 
can be neglected compared to A and to. The exchange 
of a heavy mediator with mass A may take place in the 
s-channel and/or in t-channel, as depicted in Figure 1, 
depending on the specific model. 


A. Exact formula for (av re i) 


In all generality, for 2 —> 2 processes, the matrix ele¬ 
ments depend only on two independent Mandelstam vari¬ 
ables, for example s and t, and the squared matrix ele¬ 
ment is dimensionless. After integrating over the CMF 
angle, for example, the only remaining dependence is on 
s and m. Any amplitude related to the operator ( 1 ) gives 
an integrated squared matrix element \A4 1 2 summed over 
the final spins and averaged over the initial spins that is 
a simple polynomial of the type 

ui= \M \ 2 d cos 9 = p 2 s 2 + pim 2 s + pom 11 , (2) 


with po; ■■■,P 2 depending on A and A a ,b- To get the for¬ 
mula for (av re i) in a useful form, it is convenient to define 
the reduced cross section 


<r 0 


1 1 
2 m 2 327T 


and the effective cross section 


AgA 2 to 2 


An r A 4 ’ 


( 3 ) 

( 4 ) 


which contains all the couplings. In terms of the effective 
cross section (4), and of the dimensionless variable y = 


s/(4 to 2 ), the reduced cross section Eq. (3) becomes 

cr 0 = cta (a 2 y 2 + “f ^ + Y§) > ( 5 ) 

where now a 2 , ..., a 0 are pure numbers. The total unpo¬ 
larized cross section then is 


_ 2 to 2 y/X (s,m|,m|) 

a ~ - 7= 9 = 7 = g o- (0) 

s y/\(s,m{,m%) 

We now set m\ = m 2 = to, m 3 = m .4 = 0 in 
Eq. ( 6 ) and in Eq. (B20), and change variable to y. Thus 
Eq. (B20) becomes 

2 x f°° __- 

Wei) = k2 ^ J dyy/y- lK±( 2 xffy)<To(y)- (7) 
Using the integrals of A, we find 


(av le \) = (7 a — [ 8 a 2 + 2ai + (5a 2 + 

16 


2 ai + clq ) 


KM 

K 2 {x) 


+ 3(22 


K l{x) 

K 2 {x) 


( 8 ) 


In the case m 3 = TO 4 = 0 the pure mass terms do not 
appear in the cross sections, thus a 0 = 0. Furthermore, 
we can relate a 2 and ai each other by an appropriate 
multiplicative factor, 


a\ = ka 2l (9) 

and express the cross sections as a function of a 2 only. 
The general formula ( 8 ) thus finally becomes 

Wei) = (7Aa 2 iFk(x), ( 10 ) 

with 

MX) = 16 ( 8 + 2k + (5 + + 3 ||M) 

( 11 ) 


the factored out thermal function. 

The nonrelativistic thermal average is given by the 
expansion at x 1. Using the asymptotic expansions 
Eq. (A2) we find 


( k 3 k\ 

(c7nrV r )nr = &A&2 8^* / )‘ (1^) 

In the ultrarelativistic limit, x <C 1, using the expan¬ 
sions (A3), the thermal functions behave as 3/a; 2 , thus 

{(7Vre\)ur ~ aACl2 ~2 = i^X7 3 ° 2T2 ’ 

which is the expected result for massless particles. 

The exact integration is possible because the effec¬ 
tive operator removes the momentum dependence in the 
propagators that are reduced to a multiplicative constant 
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and the assumption m 3 = m .4 = 0 allows to simplify the 
square root \J A(s, m§, m 2 ) = s in the cross section (6). 
For example, with m 3 = m .4 = m^,, equation (7) becomes 

(cru re i) = — j—r / dy y/y~~p ^/y^l K\ {2xy/y)a 0 (y, p). 

^2 \ x ) J yo 

with p = rn^/m 2 and y$ = 1 \i m > m^,, yo = p if 
m < rriip. In this case the exact integration is not possible 
but nonrelativistic expansions exist also in the case p = 1 
and p > 1 as we have shown in Ref. [27] . 


B. Applications 



In order to show the thermal behaviour of different 
interactions, we calculate the cross sections for various 
operators of the type (1), both for s and t channel anni¬ 
hilation. We list the quantity vj = A 4 /(A 2 A 2 )w and the 
resulting average Eq. (10). 

For the s-channel annihilation we find: 

1) Scalar: (xxXV’VO, (xxX^A 5 # 

n7 = 2s(s — 4?n 2 ), (asv re i) = o\lT-±{x ). (14) 

2) Pseudo-scalar: (x7 5 x)(^7 5 V’), (X7 5 X)W^) : 

^ = 2s 2 , (ap S v le i) = a A 2J 7 0 (x). (15) 

3) Chiral: (xPl,rx){^Pl,r^)- 

W = ^s(s - 2m 2 ), (<J C V rel) = CTA ^—2 (x) ■ (16) 

4) Pseudo-vector: (X7 M 75X) W’7/x75'*/X (X7 /i 75XXV’7A*V0- 


Figure 2. The thermal function (11) for the interactions and 
annihilation cross sections considered in the text. 

1) Scalar, pseudo-scalar(xx)W0> (xxX^VO. 
(X7 5 X)(V’7 5 V’) I (X7 5 X)(V’V’)■ 

w D = |s( s ~ to 2 ), ( a S,PS V rel) = CTA^J r _i(a:). (20) 

= ^s(5s — 14m 2 ), (crfpgVre l) = cta^-T 7 -!!^)- 

o ’ o 5 

( 21 ) 

2) Chiral: (xPl,rx){^Pl,r^)- 

w D = ^s(s - to 2 ), (ay?’ 4 u re i) = cr A ^ P-i(x). (22) 

= ^s(s - 4m 2 ), (cr<X’ 4 u re i) = (23) 

3) Pseudo-vector: (x7 M 75X)('07/x75^), (x7 /i 75,x)(V'7/x^)- 


8 8 

^ = 3 S ^ ~ 4m2 ^’ ( a PWre l) = &A-P- 4 (x) (17) 

5) Vector: (A H^x)^!^), (xY/xX^A/W 5 # 

8 8 

w = 3 S ^ S + 2m2 '*’ ^ avVre[S> = (18) 

6) Vector-chiral: (x7 / "1\-RX0(^7/^l,.rV’)- 

8 8 
W = gS(s - TO 2 ), (CTl/Cfrel) = CT A -.F-i (z). (19) 

The tensor interaction o tlv gives the same function as 
the vector case and is not reported. In the case of a Ma- 
jorana x clearly the vector and tensor interactions are 
absent, and the inclusion of a factor 1/2 in the opera¬ 
tor (1) cancels the factor 4 due to the presence of the 
exchange diagram of the initial identical particles. 

Now we consider some examples of t-channel annihila¬ 
tion for operators common to Dirac and Majorana DM 
annihilation: 


®D = ^ S ( 4S “ 7m2 )> ( a PV V rel) = 0"A^V_7(x). (24) 

= ^ S ( 7S “ 16TO 2 ), {apyV Te l) = CT A 16 (x ). 

(25) 


The thermal functions corresponding to the previous 
cases are shown in Figure 2 where the asymptotic be¬ 
haviours are clearly seen. In particular we note that 


Po(x) 


F-a{x) 


1 1 

(a , M{X) 


16 1 


K 2 (x 

3 1 

( K 2 {x) 

K$(x)\ 

16 1 

K KUx) 

K 2 (x)J ’ 


(26) 

(27) 


behave in the nonrelativistic limit as 

Po(x) ~ 1 + 0 (x~ 2 ), J-_ 4 (x)~^-+0( x~ 2 ). 


The function J r o(a;), which appears in the s-channel an¬ 
nihilation through a pseudoscalar interaction, is the only 
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case where the term of order 0 {x~ x ) is absent, while 
J r _ 4 (x), which appears in the scalar and axial-vector s- 
channel annihilation and in the chiral f-channel Majorana 
fermion annihilation, is the only case where the constant 
O(x 0 ) term is zero. These are the exact temperature de¬ 
pendent factors that correspond to the phenomenological 
interpolating functions proposed in Ref. [28] to model the 
s-wave and p -wave behaviour in the nonrelativistic limit. 
For all other interactions both s-wave and p -wave contri¬ 
bution are present. The function T_ 4 ( 2 ;) can be also read 
off from the formulas of Ref. [29] where the f-channel an¬ 
nihilation of Majorana fermions with the exchange of a 
scalar with chiral couplings was considered. 

We note that although we have concentrated on the 
case of fermion DM, the formula is valid for DM scalar 
and vector candidates as well, with the necessary redefi¬ 
nition of <7 A . 


III. EXPANSION OF THE CROSS SECTION IN 
POWERS OF THE RELATIVE VELOCITY 

In the general case m 3 = 7714 = m^ ^ 0 the exact in¬ 
tegration is not possible. If the relative velocity of the 
annihilating particles is small compared with the veloc¬ 
ity of light we can work directly with nonrelativistic for¬ 
mulas. The exothermic annihilation cross section in the 
nonrelativistic limit, to the lowest orders in v r , is usually 
expanded as a nr ~ a/v r + bv r , and multiplying by v r , 

<y nr v r ~ a + bv^. (28) 

Then, using Eq. (B18) and (B19), the nonrelativistic 
thermal average of Eq. (28) is 

{&nr'Ur')nr ~ U “t“ 6 . (29) 

X 

In the case of our cross sections, comparing Eq. (29) 
with Eq. (12), the coefficients are thus 1 

f k\ k 

a = a A a 2 (l + -J, b =-a A a 2 —. (30) 

We now ask, given tr(s), how to perform the expan¬ 
sion in terms of the relative velocity to find the coeffi¬ 
cients a and b that correspond to the large x expansion 


1 This result must coincide with the expansion of Ref. [27, 30]. 
With our notation the expansion is 

{&nrVr)nr ~~ ^o\y — 1 T ““ ^— <To\y — 1 T ~0'Q\y — 1^ , 

where the prime indicate derivative respect to the variable y. 
Comparison with the expansion (29) requires to identify 

a = < 7 oIb=i, b = — Cq| h= i — a o|y=l^ • 

Using Eq. (5) with ao = 0 and a 2 = kai, it is easy to verify that 
one obtains again Eq. (12). 


of the relativistic thermal average in the COF. Combin¬ 
ing equations (5), ( 6 ), (9), the general total annihilation 
cross section reads 

a 2 ( s , k \ 

,7 - <TA Ty»-4m» (irf + lj- (31) 

The correct way to proceed is to use the invariant relation 
Eq. (B5) with m\ = m 2 = m and to solve it for s as a 
function of v le \: 


,. 1 i ) ■ 02) 

V 1 - Wei/ 

This formula is valid in every frame and substituted in 
Eq. (31) gives the exact dependence of the cross sec¬ 
tion on the relativistic relative velocity, cr( 7 ! re i). Then, 
if v re i ~ v r -C 1 , we can expand the obtained expression 
to the desired order in v r and the nonrelativistic average 
taken using Eq. (B19) will coincide with the expansion 
of Eq. (10) for i> 1, that is the expansion (12). 

Equivalently, in order to find the expansion (28), we 
note that the squared roots in the annihilation cross sec¬ 
tion (31) imply that a term of order vf in s will contribute 
to the order d 2 in a. Thus we need to expand s, formula 
(32), at least to order vf, 

s ~ 4to 2 + to 2 ?; 2 + ^to 2 vf. (33) 

Substituting Eq. (33) in Eq. (31) and performing the ex¬ 
pansion in powers of v r it easy to find 

f k k \ 

a nr v r ~ a\a 2 ( 1 + - - — J , (34) 



in agreement with (30). 

In the case of coannihilations [31], for example when 
a DM particles scatter off another particle with different 
mass, the Mandelstam invariant takes the form 


s = (toi — m 2) 2 + 2 toiTO 2 1 + 


V 1 - v l 


rel, 


with the expansion 


s ~ (mi — m 2) 2 + mim 2 vl + -mim 2 vf. 


(35) 


(36) 


This procedure gives the correct expansion in the COF 
where the velocities Vi 2 of the colliding particles are 
specified in this frame. Clearly, the same expansion with 
the same coefficients is obtained in the LF and in the 
CMF. 


IV. THE PROBLEMS WITH THE M0LLER 
VELOCITY 

The simple outlined procedure has not been recognized 
in DM literature where, incorrectly, the Mqllcr velocity 
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v. Eq. (B22), instead of D re i is considered. As reminded 
in B, v is a non-invariant, non-physical velocity. The 
expression of v in terms of s is thus different in different 
frames and the expansion of a takes different values in 
different frames. 

Before discussing the problems with the Mpller velocity 
we note that if we take the limit m/ —> 0 in the analogous 
expansions published many papers [4-11], we do not 
reproduce the expansion (34) . The reason is that in these 
papers the expansion of s is truncated to the lowest order 
in v%, 

s ~ 4m 2 + to 2 d 2 . (37) 

If we substitute this in Eq. (31) and expand, we find 
( k 12 + k 2 \ 

CTnrVr ~ 0^02 I 1 + - H- I > (38) 


with an incorrect coefficient b. Clearly the same wrong 
result is obtained truncating (33) to order v%, whatever 
the frame in which v r is specified, CMF, LF or COF. 

We now go back to the Mpller velocity (B22). Evalu¬ 
ated in the CMF taking m\ = m 2 = to reads 


D* 


— 7 =\/s* - 4m 2 . 
v®* 


(39) 


We indicate the quantities evaluated in the CMF with a 
By inverting Eq. (39) we find 


s* — 



(40) 


This relation is different from (32) and is often incorrectly 
identified as the relation between s and the relative ve¬ 
locity in the CMF, see for example [31], [10]. In facts, 
the expansion to order reads 


s* ~ 4m 2 + to 2 !) 2 * + -^-!>)f*. (41) 

When used in (31), it gives the following nonreltivistic 
expansion of the cross section 


O’ nr'Or ,* 


CTA02 1 + - + 


(42) 


which is different from the correct expansion (34). 

Other authors, [26] and [20-23], perform the expansion 
with the Mpller velocity evaluated in the rest frame of one 
particle. Indicating with the quantities in this frame, 
Eq. (B22) becomes 


Vi 


y/sIVst - 4to 2 
Si — 2 to 2 ’ 


and by inverting Eq. (43) we obtain 


Si = 2 to 2 



(43) 


This expression is formally identical to Eq. (32), thus 
when Vi ~ v r> e and se is expanded up to the order id t 
we obtain the expansion o nr v r / which formally coincides 
with Eq. (34), with v rt i in place of v r . 

It should be clear that this is just a mathematical 
coincidence due to the fact that v reduces to !) re i only 
when one of the two velocities iq .2 is zero as it is evident 
from the definitions Eq. (B2) and Eq. (B22). In other 
words, the expansion found in Refs. [20-23] are correct 
because the authors have implicitly used the relative ve¬ 
locity, Eq. (B5) and (33). 

We thus emphasize some common statements found in 
DM literature and why they do not subsist: 

1) In the relativistic Boltzmann equation the v in av is 
v and (av) must be calculated in the LF frame. 

This is not true, as shown in details in Ref. [27] and in 
B. Using ii re i and recognizing the nonphysical nature of 
v. one works always with invariant quantities and the 
consistency of the relativistic and nonrelativistic formu¬ 
las and expansions is obtained in the comoving frame 
without any further specification of the kinematics. The 
LF, also called Mpller frame in Ref. [23], cannot be a 
privileged frame for the relic abundances calculation also 
because for massless particles the rest frame does not ex¬ 
ist. 

2) The M 0 ller velocity coincides with relative velocity in 
a frame where the velocities are collinear. 

This not true because, for example, in the CMF where 
the particles have velocities d* , the Mpller velocity is 2v * 
while the relative velocity is 2v*/(l + d 2 ). Note that the 
true relative velocity is never superluminal. 


V. IMPACT ON THE RELIC ABUNDANCE 

Only in the case k = —4 the incorrect expansions (38) 
and (42) coincide, incidentally, with the expansion (34). 
While the lowest order coefficient a turns out to be al¬ 
ways the same, the coefficient b is different in any other 
case. To illustrate the impact of b on the value of the relic 
abundance we consider the case of the s-channel annihila¬ 
tion with vector interaction, Eq. (18), and the s-channel 
annihilation with a pseudoscalar exchange, Eq. (15). In 
the first case k = 2, 02 = 8/3, and the correct coefficients 
a and b are 

ay = 4<t a , b v = -^, (45) 

while the incorrect coefficient b in (38) and (42) is 
7 2 

bv1 = e CT A, b V2 = — <t a . (46) 

6 3 

In the second case, k = 0 and <12 = 2, thus 

aps = 2cta, bps = 0, (47) 

and the wrong b coefficients are 

bps-L = ^ a , bps 2 = -a\. 


(44) 


(48) 
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We calculate the relic abundance following the exact 
theory of freeze out presented in Ref. [32] . We briefly re¬ 
call the main points. Let Yq = 45/( 47 r 4 )(g x / g s )x 2 K 2 (x) 
be the initial equilibrium abundance (number density 
over the entropy density), with g x = 2 for spin 1/2 
fermions and g s the relativistic degrees of freedom as¬ 
sociated with the entropy density. The function Y\ (x) 
that gives the abundance up to the point x* where 
Yi(x) — Yq(x) is maximal is 


Fi(x) 

5{x) 


(1 + S(x))Y 0 (x), 



x 2 1 dY 0 
C(av le \)Y 0 Y 0 dx 


- 1 , 


(49) 

(50) 


with x* given by the condition 


1.12 


1.10 


1.08 


1.06 


1.04 


1.02 


1.00 



100 200 300 400 500 600 700 800 9001000 100 200 300 400 500 600 700 800 9001000 

m (GeV) m (GeV) 


_1 dY 0 (x) _ 1 dS(x) f51) 

Y 0 (x) dx S(x) dx *' 

The abundance at x > x* is found by integrating numer¬ 
ically the usual equation 

rlV C 

— = ^(av iel )(Y 2 -Y 2 ), (52) 

with the initial condition (x*, V(x*) = Yi(x*)). The fac¬ 
tor C is defined by C = Mpm xy fgl , where Mp is the 

Plank mass and yfgZ = g s /yfgp{ 1 + T/3d(\ng s )/dT) ac¬ 
counts for the temperature dependence of the relativistic 
degrees of freedom associated with the energy density, g pi 
and g s [26, 30]. For WIMP masses larger than 10 GeV we 
can neglect the temperature dependence of the degrees of 
freedom [33, 34] and take g s = g p = g = 100, yd/f = yfg- 
In solving numerically (52) and (51) with the exposed 
method, we use the exact formula for (av le i), Eq. (10). 

We compare the previous numerical solution with the 
one obtained using the nonrelativistic freeze out approx¬ 
imation (FOA) that is commonly employed in literature. 
The FOA consists in integrating equation (52) with an 
initial condition ( Xf , Y(xf)) such that the equilibrium 
term proportional F 0 2 can be neglected. We choose the 
freeze out point at the point X 2 where Y{x 2 ) ~ Yi(x 2 ) = 
2 Fo(x’ 2 ). As shown in Ref. [32], Yi(x) well approximates 
the true abundance also in the interval x* < x < X 2 . X 2 
is the optimal point for the FOA and corresponds to the 
temperature where the extent of the inverse creation re¬ 
action iptjj —> XX is maximal. The solution in the freeze 
out approximation is then 


v _ _ 2Yq{x 2 ) _ /rox 

FOA l + 2Y 0 (x 2 )£(a + 3±y 

The freeze out point x 2 is given by the condition — ^ 

= 3^(f7n re i)Fo, which, in terms of the method of Ref. [35] 
corresponds to c(c + 2) = 3, that is c = 1. Using the 
nonrelativistic form of Yq, 


45 9x 


4tt 4 g s 



(54) 


Figure 3. Ratio of the relic abundance obtained by solving nu¬ 
merically equation (52) over the value given by the freeze out 
approximation, for the pseudoscalar and vector interactions. 
In the bottom blue curves for the FOA the correct coefficients 
(45) and (47) are used. The red and the black curves show 
the effect of the wrong coefficients (46) and (48), respectively. 


x 2 is given by the root of 

3C + 6 —^ f-x^e-* = 1. (55) 

Calling a = 3aC^/n/2, an accurate analytical approxi¬ 
mate solution of Eq. (55) is given by 

1 /? 7 

X 2 = lncr-ln(lna) + ln(l -|-(InaW 1 ). (56) 

2 a 

The relic abundance normalized over the critical den¬ 
sity is flh 2 = 2.755 x 10 8 (m/GeV)Y( oo ) for a Majorana 
fermion and two times that quantity for a Dirac fermion 
with the same density of antiparticles. We now compare 
the exact relic abundance Qh 2 with the value ( tth 2 )poA 
furnished by the nonrelativistic FOA calculated using the 
correct and the wrong expansions. We take the couplings 
A a ,b = 1 for illustrative purposes and two values of the 
cut off scale, A = 1,10 TeV. The value of the freeze out 
points x* and X 2 varies roughly between 18 and 30 in the 
parameter space with m < A where the effective treat¬ 
ment is supposed to be valid. The ratio Q.h 2 /{VUi 2 )foa 
is shown in Figure 3 as a function of the DM mass for 
the chosen examples. The bottom blue curves show that 
the FOA with the correct coefficients (45) and (47) un¬ 
derestimates the numerical value by less than 2 %, and 
that in most part of the parameter space the error is at 
the level of 1% or less. This a test of goodness for our 
FOA, and confirms what shown in Ref. [32]. The red 
and the black curves show the effect of the wrong coeffi¬ 
cients (46) and (48), respectively. The wrong expansions 
underestimate the relic abundance by a factor between 
3% and 12% for both interactions for masses larger than 
10 GeV as shown in the plot. The behaviour is similar 
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for the other interactions not shown in figure. The er¬ 
ror becomes even larger at smaller masses and we have 
verified that using for example c = 1/2 and other values 
we get even worst approximations. Clearly this kind of 
error nowadays is not compatible with the precision with 
which the experimental value is known. 
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Appendix A: Integrals and expansions 


Equation (7) can be written as 

( ™ rel) = CTA i^) ( azA2 + T Al + S A °) ■ (A1) 


The integrals are evaluated with methods similar to those 
described in Ref. [27] in terms of Bessel functions of the 
second kind: 


Ao — 

Ai = 

A 2 = 


/ °° _ i 

dy\/y — IK^Xyfij) = 

f°° , ,-- T ., „ 1 K 2 (x) + K%(x) 

J^ dyyjy- lyKi(2xy/y) = — --—=-. 

/ oo _ 

dy\Jy- ly 2 K 1 (2x^/y) 

±-±[5K 2 (x) + 8 Kl(x) + 3 K 2 (x)}. 


The expansions at x 1 are 


KM 

K%(x) 


l--+0 (x~ 2 ), 

X 


Kj(x) 

Ki(x) 


1 + -+0(x~ 2 ), 

X 

(A2) 


while for x <C 1 are 


km 

I< 2 {x) 


x 2 

T 


0 (x 3 ), 


Kj( x) 
K 2 {x) 


1 

-+0(x 2 ). (A3) 

x z 


Appendix B: Invariant formulation using v re i 

In this Appendix we remind, based on the results of 
Ref. [27], the main points about the relation between the 
relative velocity, the Mpller velocity, flux and thermal 
average which are used in the main text. 


1. Invariant relative velocity 

The relativistic relative velocity that generalizes the 
nonrelativistic relative velocity 

v r = \vi-v 2 \, (Bl) 


is given by 


^rel 


y/{v i - v 2 y - 

1 _ Vi-V 2 


(B2) 


We have explicitly written the dependence on the velocity 
of light c to make manifest that v le \ coincide with v r in 
the nonrelativistic limit because the scalar and vector 
products are of order {v/c) 2 . In the following we go back 
to natural units. 

The relative velocity zy e i can be written using the Man¬ 
delstam invariant s = (pi +p 2 ) 2 , where pi j2 are the four- 
momenta, and A, the Mandelstam triangular function, 

\{s,m\,m\) = [s — (mi + m 2 ) 2 ][s — (mi — m 2 ) 2 ], 

(B3) 

in a generic frame, 


(B4) 


(B5) 

showing its invariant nature. 


_ V (Pi ■ P2 ) 2 - mlml 
^rel — 

Pi ‘ P2 

= VA(s,mj,m^) 
s — (m 2 + m\)' 


2. Flux factor 

Given two bunches of particles with number densities 
ni >2 and velocities iq j2 in a generic inertial frame, in non¬ 
relativistic physics the flux is F nr = n\n 2 V r . To obtain 
the relativistic invariant flux that reduces to F nr in the 
nonrelativistic limit, the easiest way is to consider the 
4-currents Ji = (ni,mui), thus 


F = (Jr • J 2 )w re i = nin 2 (l - v 1 ■ v 2 )v iel . (B 6 ) 


Note that the factor (1 — i?i • v 2 ) that guarantees the 
Lorentz invariance of the product of the number densities 
can also be written as 


1 — Vi ■ v 2 = 


7r 

7172 


Pi ■P2 

e x e 2 ' 


(B7) 


where y r = l/i/l — v 2 el is the Lorentz factor associated 
with i> re i and 7 j the Lorentz factors associated with 17 . 

If the element of Lorentz invariant phase space is de¬ 
fined as usual 


- _ d 3 p t 
Vl (2tt) 3 2 Ei ’ 


(B 8 ) 























and one particle states for bosons and fermions are nor¬ 
malized to 2 Ei such that the density per unit volume is 
2Ei, then, using (B7), the flux (B 6 ) simplifies to 

F = 4(pi -p 2 )v le i. (B9) 

Substituting the expression of u re i in the momentum rep¬ 
resentation, formula (B4), in Eq. (B9), the scalar product 
Pi • P 2 cancels out and the standard explicit form is re¬ 
covered 

F = 4^/(pi • P 2) 2 — m\m\. (BIO) 

3. Cross section and collision integral 


When the species 1 and 2 are the same, it takes the usual 
form (av re \)(riQ — n 2 ) with the factor 1/2 cancelled by sto¬ 
ichiometric coefficient appearing in the left-hand side of 
the complete kinetic equation, see for example Ref. [32]. 

4. Averaged thermal rate 

In Eq. (B13) the general definition of relativistic ther¬ 
mal averaged rate is 

r 2 d 3 p 

We l) = / n^ 1 •P 2 )/o,p(Pl)/o,p(P 2 ) <7 ^rel- 

J i= 1 1 

(B14) 


The integrated collision term of the Boltzmann equa¬ 
tion, neglecting quantum effects, can be written as, 


In the case of the relativistic Maxwell-Boltzmann-Juttner 
statistics, the momentum distribution is 


/ n d Pi [/ 3 / 4 W3W,2) Wi/2W(1,2|3,4)], 

^ i= 1 

where W(ij\kl) = (2t r) 4 <5 4 (Py - Pki)J2 Si , Sf l-WW 2 > 
and fi is the phase space distribution. 

Using the unitary condition f dp 3 dp4W(3, 4| 1, 2) = 
/ dp 3 dp 4 ,W(l, 2|3,4) to write the collision integral only 
in terms of the annihilation rate 

1 f ^ 

YW— J IT d & (/a/4 - hh)W(l, 2|3,4), (Bll) 


we keep out a statistical factor accounting for the possi¬ 
bility of identical particles. 

By definition, the invariant cross section, using the flux 
in the form (B9), is 


1 


4(pi -p2Hel 


dp 3 dp 4 


WW2|3,4) 

9i92 


(B12) 


being g l = ( 2 s^ + 1 ) the spin degrees of freedom. 

Assuming as usual that the annihilation products are 
described by the equilibrium phase space distribution at 
zero chemical potential f oi , we have f 3 f 4 = /o, 3 /o ,4 = 
/o,i/o, 2 ; the last equality following from energy conser¬ 
vation. Hence 


/o,p(p) 


1 c -\/p 2 +m 2 /T 

AirirPT K2 {x) 


(B15) 


and as shown in Ref. [27] , the six-dimensional integral on 
the right-hand side of Eq. (B14) reduces to 


(crurei) = / dv le \V(v le ])av le i, 


(B16) 


where the probability distribution of i> re i, for example for 
mi = m 2 = m, is 


P(u re i) 


x 7 r 3 (7 2 -l) 

V2Kl(x) VtT+T 


+ 1 ). 

(B17) 


This is completely analogous to the nonrelativistic case 
where the probability distribution of v r , for mi = m 2 = 
m, is 


P(v r ) = \ -x 3/2 v 2 e~ x ^, 


(B18) 


and the thermal average reads 


/»00 

(^nr^r)nr = I dv r P(v r ') 0 ’ nr V r . (B19) 

Jo 


J LI • ^)(/0,l/0,2 - /l/2)^rel. 

The equilibrium phase-space distribution fog is related 
to the number density no and to the momentum distri¬ 
bution /o,p(p) by /(27r) 3 fog = n 0 gf 0 , P (p)■ Assuming 
further that the non-equilibrium phase-space function at 
finite chemical potential /j remains proportional to the 
equilibrium momentum distribution by a factor given by 
the non-equilibrium number density ni, gi/(2n) 3 fi = 
nifo, P (p), we obtain 

(n 0 ,mo ,2 - nin 2 ){crv ie i). (B13) 

1 + C>12 


Given the total annihilation cross section a the product 
av re \ will reduce to the nonrelativistic limit a nr v r and 
(tTUrei) to ( CT nr v r )nr in the COF when v re i ~ v r -C 1. 
Expressing Eq. (B16) in terms of s using (B5), we obtain 
the usual integral [27, 36] useful for practical calculation 


(av le i) = 


srrLmWW 


f°° A(s, ?nf,m 2 ) ,^7*, 

/ ds - 7 =- Ki(—)a, 

Im 2 V s 4 


(B20) 


with Xi = rrii/T and M = (mi + m 2 ). 

We have recently become aware of the paper [37] 
where, probably for the first time, the thermal average of 
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relativistic rates was discussed and it was realized that 
with the relativistic Maxwell-Boltzmann statistics for¬ 
mula (B14) reduces to a single integral over the distri¬ 
bution over the relative momentum. With some algebra 
and change of variables it is easy to verify that for exam¬ 
ple Eqs. (lib) and (12a) of [37] coincide with Eqs. (29) 
and (37) of Ref. [27]. In Ref. [37] the cases of collisions of 
two massive particles, two massless particles and a mas¬ 
sive with a massless particles are treated separately as if 
different definitions of flux and cross sections were neces¬ 
sary in each case. Clearly this distinction is unnecessary 
for the formulation we have given is completely general 
and valid in any case. We finally note that an integral 
formula similar to (B20) was also given in Ref. [29]. 


5. No need for the Mpller velocity. 

By noting that in Eq. (B6) the factor (1 — iq • v 2 ) can 
cancel the same factor in the denominator of v re i, the 
invariant flux can also be written in the form 

F = mn 2 V^i ^ v z) 1 2 3 4 5 6 7 - ( v i x v 2) 2 - (B21) 

In the textbook by Landau and Lifschits [38] this form 
is attributed to Pauli without giving any reference, while 
its origin is more generally attributed to Mpller [39] . 

It is interesting to look at original paper by Mqllcr [39] . 
With our notation, he wants to prove that the flux given 
(B21) is invariant. In order to do that he shows that this 
can be written as a product of two in variant quantities: 
the ratio and the quantity B = 1 J(pi ■ p 2 ) 2 — m\ 
and there he stops. 


The flux factor written in the form (B21) has the same 
structure of thee nonrelativistic expression n\n 2 v r . Prob¬ 
ably for this reason it has been later introduced in the 
literature the notion of Mpller velocity 


v = \f{v 1 - v 2 ) 2 - (v 1 x v 2 ) 2 


V(pi T2) 2 ~ mjml 
EiE 2 

(B22) 


It is worth to stress that neither Mpller nor Landau and 
Lifschits attribute any particular meaning to Eq. (B22) 
and do not define it as a particular velocity, even less as 
relative velocity. Clearly v is nothing but the numerator 
of the formula defining u re i because v = (1 — v-± • v 2 )v le \, 
where the factor (1 — V\ • v 2 ) comes from the definition of 
the invariant flux (B6). Already this fact indicates that 
v is not a fundamental physical quantity and overall, it 
is not the relative velocity, nor when the velocities are 
collinear. 

On the contrary, in DM literature and in textbooks, 
when defining the flux factor for the relativistic invariant 
cross section, it is incorrectly asserted that in a frame 
where the velocities are collinear the quantity \vi — v 2 \ is 
the relative velocity, while in a generic frame is given by 
(B22). The form (B21) of the flux is a simple consequence 
of the fundamental quantities (B2) and (B6), there is no 
new physics or concept in it. For these reasons, and for 
its noninvariant and nonphysical nature, v should not be 
used. 
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